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Abstract: 

We show that the conductance (G) and transition voltage (V,) spectroscopy data for standard vacuum nanogaps reported in the 
preceding paper by Sotthewes et al cannot be understood within the framework of the existing theories/models whatever realistic 
effects are incorporates. However, if we include an additional (“ghost”) current, the trends of the dependencies of G and V) on 
the nanogap size d can be explained. The ghost current is very small. Therefore, effects related to the ghost current can only be 
revealed at larger c/’s, where the ghost current becomes comparable to or overcomes the tunneling current. Although we are not 
able to unravel the origin of this ghost contribution, we can and do refer to experimental G- and VJ-data reported for molecular 
junctions, which exhibit similarities to the presently considered vacuum nanojunctions. Analyzing notable difficulties related to 
the regime assigned as transport via hopping, we speculate that a ghost contribution may also exist in molecular junctions. 
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1 Introduction 

A long standing challenge in the field of molecular electron¬ 
ics is whether the measured currents are really mediated by 
the active molecules intentionally embedded in the nanode¬ 
vice under consideration or are (also) significantly affected, 
e.g., by impurities or defects.^ From this perspective, the 
transport through vacuum nanodevices seems to be privileged, 
because it is hard to imagine how could the measured current 
I = Imeasured be plagued by undesired “ghost” channels giving 
rise to a current Ighost overimposed on the current Inmnei due to 
the tunneling across the vacuum nanogap. However strange 
it might be, after ruling out all causes we could conceive, 
we have to conclude that the challenging experimental results 
on the transition voltage V) spectroscopy (TVS)^ reported in 
the preceding paper— cannot be explained without assuming a 
“ghost” contribution. The origin of such a transport channel 
contributing to the transport in standard vacuo remains the big 
mystery, which the present work cannot unravel. Still, what 
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we can presently do is to point out to some similarities with 
results reported for molecular junctions^”— and suggest that a 
similar contribution can also be present there. 

2 Theoretical details 

The tunneling current through a vacuum nanogap can be de¬ 
termined by exact numerical solution of the Schrodinger equa¬ 
tion by using a tunneling barrier 0^ that includes contribution 
from the applied bias V, electrodes’ (STM substrate s and tip 
t) work functions <I>sand image charges (j)im^ 

Mr,V) = (<I>f+<I>,)/2 + 0;„,(z)+yfc(z), (1) 

Vh{z) = {<^s-<^t-eV)z/d. (2) 

The “bare” rectangular barrier due the electrodes’ work func¬ 
tion is tilted by the applied bias V and reduced by (classical) 
image charges For two planar electrodes, the exact 

expression of 0,„,— can be used to numerically find the exact 
transmission coefficien t 

In the discussion that follows we will assume that, besides 
the current Inmnei due to electron tunneling through a metal- 
vacuum-metal junction in a scanning tunneling microscope 
(STM) setup, the measured current also comprises a “ghost” 
current 4/,os, (c/[B 

I — Imeasured ~ Itwmel Ighost- ( 3 ) 
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Fig. 1 Schematic representation of the transport scenario proposed 
in the present paper. 


For simplicity, we will assume that the latter represents an 
Ohmic contribution, that is, 

hho.t = ^. (4) 

Here V is the applied bias between the STM tip and substrate, 
which are spatially separated by a vacuum nanogap of size d, 
and p stands for the ghost (average) resistivity. 

Below, we will focus our attention on the Ohmic conduc¬ 
tance G and the transition voltage V). V) is defined as the bias at 
the minimum of the Fowler-Nordheim quantity ln(//y^)— , or, 
equivalently, by the point where the differential conductance 
is two times larger than the pseudo-Ohmic conductance^. 
The interplay between the tunneling and ghost contributions 
of eqn Q will be characterized by the dimensionless quantity 
p related to the ghost resistivity p 

p~ 1.618 (5) 


3 Results 

3.1 Attempting to incorporate realistic effects into the 
tunneling barrier 

The results for V) presented in Fig. 1 and 5 of ref. 0, showing 
that calculations based on Simmons-type approximation, simi¬ 
lar to those done earlie r cannot explain the tf-dependence 
of the transition voltage V) found in experiment {cf. Fig. 4A 
and 4B of ref. 0) may not be surprising; shortcomings of the 
Simmons approach to transition voltage spectroscopy are am¬ 
ply documented '*^4^ . 

In the first part of the present analysis, we will show that 
the tunneling transport model cannot explain the experimen¬ 
tal V)-data even when the numerical exact solution of the 


Schrodinger equation is employed (this is the case for all re¬ 
sults presented here), and even attempting to make the barrier 
model more realistic. This is shown by the results presented in 
Fig. S1 of the ESIf, which emerged from the attempt to make 
the tunneling barrier as realistic as possible. (Throughout, la¬ 
bels S refer to the ESIj.) 

Eig.|2^ shows the effect of the difference between the STM- 
tip and substrate work functions (<t>, —Ovr = 4.55 eV versus 
<I>i <!>/>, =5.65 eV or <I>, 4 u =5.2 eV, respectively, cf. eqn O- 
This may be important, because it yields an extra Volta-type 
field of a strength |<I>i — <I>f|/(ec/) ^ 1 V/nm comparable to 
the applied field V As seen in Eig. |2^, the asymme¬ 
try semi-quantitatively explains the small difference 

Vt+ f IV-I between positive and negative biases for Ft sub¬ 
strates. As noted—, the Shockley surface state, not included 
in the present calculations, represents a plausible source for 
the considerably larger asymmetry \Vt-\ for Au (111) 

substrates. Eig. Sla demonstrates that image effects, which 
may be important at smaller sizes^, are not essential for 
the values used in experiments—. The same figure indi¬ 
cates that the potential profile has no significant impact on 
Vt (cf. Eig. Sla). Instead of well defined d- and O-values, it 
might be more realistic to consider statistical distributions of 
these quantities. These distributions may better model STM- 
tips without a well defined facet and polycrystalline Pt sub¬ 
strates. But, as visible in Eig. Sib and Sic, these effects do 
not improve the agreement with experiment. 

Eig. |2]3 depicts results for the Ohmic conductance G, a 
quantity that also plays an important part in the subsequent 
analysis. Whatever the effects considered (cf. Eig S2), it does 
not notably deviate from the well-known exponential decay. 
This is also at variance with the experimental G-data, which 
evolve to a plateau at larger ^f-values (cf. Eigure 4C of ref 0 ). 

3.2 Impact of a ghost current 

Let us now assume that a ghost current [eqn (|4|] exists, and ex¬ 
amine its impact on G and V). Results for G and V) at various 
values of the parameter p are presented in Eig.[^ and Eig.j^t. 
The trend observed in these figures is compatible with the ex¬ 
perimental behavior. 

In the plateau regime (d ~ 1.12nm), experiments found 
Gpt ~ 5 X 10^® S for platinum substrates and Gau — 5 x 
lO^'^’S for gold substrates (cf. Eigure 4C of ref. 0 ). A reli¬ 
able estimate for the transverse area vacuum nanogap depends, 
amongst others, on the actual tip shape and cannot be given—. 
So, let us assume a tip radius r = 1; 10; lOOnm. Inserting the 
above values in eqn (|5]), we get p ~ 10^; 10^; 10® for Au and 
p ~ 10"^; 10®; 10^ for Pt. These estimates fall in the p-range 
employed in Fig. [3] wherein the increase of V) with 1 /d pre¬ 
dicted by the conventional theory switches to the decrease ob¬ 
served in experiments^. So, a ghost current may be at the ori- 
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Fig. 2 (a) Transition voltage Vt and (b) Ohmic conductance as a 
function of the nanogap size d for STM-tip of tungsten 
(<f>iy = 4.55 eV) and substrates of Pt (<Fpf = 5.65 eV) and Au 

= 5.2 eV). Image charge effects are included using the exact 
classical interactio n which is cutoff close to electrodes using 

the procedure described elsewhereJi. 


gin of both Ohmic conductances saturating at larger d's and 
transition voltages that decrease with 1/d. We have checked 
that this trend is not affected if the realistic effects discussed 
in the preceding subsection and the ESIj are included in cal¬ 
culations. 

As a supplementary test related to the most unexpected ex¬ 
perimental finding of ref. 0 (namely, a Vt decreasing with 
1 /d), we present further results in Fig. |4] The black squares 
shown there are experimental data Vt vs. 1 /d. The other sym¬ 
bols depicted in Fig.|4]correspond to Vt obtained by subtracted 
from the experimental I — V curves a (ghost) current propor¬ 
tional to V [cf. eqn (IHi], which has been gradually increased 
from cyan-magenta-blue-green to red circles (/g/,o 5 t-values are 
indicated in the legend). In spite of certain scattering in the 
data, the trend is pretty clear: the V) data switch from a V) 


decreasing with 1 /c/ to a V) increasing with l/d. The differ¬ 
ences between the red circles and the red line (linear fit of 
the red circles) of Fig. 0] are within the experimental errors. 
We assign the last curve (red circles) as a situation where the 
ghost current has been eliminated. So, what remained is the 
genuine tunneling current, and the corresponding curve (red 
points) Vt vs. 1 /d behaves like the curves found in Fig. |2^ 
within the conventional theory, wherein merely the tunneling 
current has been accounted for. The red circles correspond 
to a ghost resistance Rghost ~ 10 GO. Assuming again tip radii 
r = 1; 10; 100nm, we find p 0.4 x 10®;0.4 x 10*;0.4 x lO***, 
which is compatible with the range covered in Fig. [3 

To end this section, we note that neither the d- nor the In¬ 
dependence of the ghost current assumed in eqn (|4|i is critical 
for the behavior depicted in Fig. [3 To check the impact of 
the 1 -dependence of eqn (|4]i, we performed calculations by 
using a -independent ghost resistance and found no signif¬ 
icant modification. Fikewise, to check the impact of the In¬ 
dependence, instead of the contribution of eqn (|4]i, we added 
a small constant term T to the transmission across the tunnel¬ 
ing barrier [i.e., using 3^{E^-,V) £^{E^-,V) -b T in eqn (SI)]. 

Results obtained in this way, which are presented in Fig. S3, 
do not qualitatively differ from those of Fig. [3 Notice that 
this procedure — which amounts to consider a “ghost” trans¬ 
mission channel, a pseudo-diffusion or hopping — yields via 
eqn (SI) a nonlinear ghost current. To summarize the results 
based on further functional dependencies studied, we can state 
that the qualitative behavior depicted in Fig. [3holds provided 
that the d- and V-dependencies of the ghost current are much 
weaker than the corresponding (exponential) dependencies of 
the tunneling current. 

4 Discussion 

We have shown above that the experimental data for electron 
transport in STM vacuum junctions can be rationalized by pos¬ 
tulating the existence of a ghost current. What could be the 
origin of the ghost current? A imperfect electric insulation 
cannot come in question—. The vacuum in experiments is 
very good (< 3 x 10^*^ mbar); therefore it is very unlikely that 
a conduction channel through the vacuum junction is respon¬ 
sible for the ghost currents needed in Fig. [3 and Fig. 01 To rule 
out a ghost current flowing through the electronic circuit, the 
experimentalists have checked that for large tip-sample sepa¬ 
ration (^ 2nm) the measured current is negligible-!^. A small 
DC offset voltage of the I — V converter gives an unavoidable, 
small contribution (^1—2 pA) to the measured current, which 
should be independent on the bias; however, any residual off¬ 
set current has been subtracted before processing the data—. 
Therefore, this possibility should also be ruled out. 

So, we are unable to unravel the origin of the ghost cur¬ 
rent postulated above (even, see below) in a standard vacuum 
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nanogap. Still, we believe that this issue may be related to 
pieces of experimental work on molecular junctions of two 
independent groups. Some time ago, Frisbie’s group first re¬ 
ported a crossover between two different transport regimes in 
molecular junctions employing the conducting probe atomic 
force microscopy (CP-AFM) setup^. For shorter molecules, 
the resistance R was found temperature (T) independent and 
increasing exponentially with increasing length (d), in con¬ 
trast to longer molecules characterized by an Arrhenius T- 
dependence R exp{Ea/kBT) and a linear c/-dependence. 
Extracting information on single-molecules using transport 
data on CP-AFM junctions^, comprising many (~ 100) 
molecules, might be not straightforward. A subsequent work 
of Tao’s group-ion single-molecule STM-junctions confirmed 
the crossover between two transport regimes with the same 
signatures noted above^. Fully in light with current under¬ 
standing of charge transport, these d- and T -dependencies 
have been taken as clear evidence of charge transport via 
tunneling in shorter molecules and via hopping in longer 
molecules^"—. 

While by no means intending to challenge that, with in¬ 
creasing d, tunneling becomes less effective and hopping in¬ 
creasingly important, for reasons exposed below we believe 
that this is not the whole issue related to the transport data of 
refs. iS and 0 Particularly problematic to us appear the 
very large values of the activation energies extracted from ex¬ 
periments, an aspect on which ref. @ already drew attention. 
As discussed in the ESIf, ref. 0 drastically overestimated the 
reorganization energies X, because the contributions of all in¬ 
tramolecular vibrations were included in calculations; in real¬ 
ity, only very low frequency modes, which can be thermally 
activated, can contributed^. Even though, the (over)estimated 
A-values remain several times smaller than those needed to ex¬ 
plain the large experimental activation energies Ea- The mea¬ 
sured values Ea — 0.54 — 0.62 eV— and Ea — 0.58 — 0.58 eV— 
would require enormous reorganization energies il > 2eV 
[cf. eqn (S5)] hardly compatible with known (including au¬ 
thor’s) results of quantum chemical calculations. 

The (f-dependence of the resistance R of the aforementioned 
experimental works on molecular junctions^“-d is similar to the 
presently considered vacuum nanojunctions (compare Eigure 
2A of ref. 0 Eigure 6b of ref. 0 and Eigure 3 of ref. 0 with 
Eigure 4C of ref. 0 and the present Eig. [^j). This represents 
one important fact for the present analysiSjbut not the only 
one. In addition to results for R vs. d, ref. Is] reported results 
on Vt for oligophenyleneimines (OPIs) of various sizes. Inter¬ 
estingly, the trend visible in Pig. |5l wherein these results are 
depicted, resembles the trend of Pig. [^. So, it is tempting 
to interpret the similarity between the ^-dependence of Pig. |5] 
and Pig. as suggesting that a ghost current may also exist 
in molecular junctions. 


5 Conclusion 

TVS was initially proposed as a tool of extracting the energy 
offset of the dominant frontier orbital relative to electrode’s 
Permi energy in molecular junctions^, and a series of stud- 
demonstrates that it holds its promise. The fact 
that Vf can be used to quantify a significant nonlinear trans- 
port-i^ makes TVS a method of interest even beyond molec¬ 
ular electronics, and ref. 0 as well as the present paper have 
illustrate its utility for vacuum nanojunctions. 

We have shown that the conductance and TVS data for a 
standard vacuum nanogap reported in ref. 0 cannot be un¬ 
derstood within the framework of the existing theories/models 
whatever realistic effects are incorporates. However, if we in¬ 
clude an additional (“ghost”) current, the experimental trends 
G ws. d and V) vs. d exhibited by the Ohmic conductance 
and transition voltage can be explained. The ghost current 
is very small; it could not have played a significant role in 
ref.— because at the very small nanogap sizes employed there 
(d < 0.4 nm—) the tunneling current is too large, and it masks 
a possible ghost contribution. Effects related to the ghost 
current can only be revealed at larger nanogap sizes d (like 
those of ref.—), where the ghost current overcomes the tunnel¬ 
ing current. This is the reason why we believe that incorpo¬ 
rating electrodes’ atomistic structure {e.g., within approaches 
based on the density functional theor y will not change 
the conclusions of the present study, which is mostly based 
on the classical barrier picture for charge transport via tunnel¬ 
ing in vacuo. Such details can be quite relevant for TVS in 
other cases: e.g., for vacuum nanogaps of sizes smaller than 
those {d > 0.8 nm) presently considered—, or molecular 
tions^, where through-space processes are negligible. 

The analysis of existing studies on size-dependence of G 
and Vt in molecular junctions have indicated that a ghost cur¬ 
rent may also be present in those nanosystems, which could 
have a significant role in properly assessing and characteriz¬ 
ing the transport mechanism (e.g. hopping). Erom this per¬ 
spective, investigating the impact of temperature on G in the 
non-exponential c/-dependent regime is a possible topic of 
further experimental work. Another subject of experimen¬ 
tal investigation could be the shot noise. Puzzling, although 
shot noise is a concept put forward in the context of vacuum 
tubes^, we are not aware of shot noise ,5^,/, studies in vac¬ 
uum nanogaps. Such studies may provide supplementary in¬ 
formation. Both G and V) sample information contained in the 
I — V -characteristic (expressed via the first power of the trans¬ 
mission coefficient ^iT), while also comprises information 
on the second power in addition, at relatively large biases 
V ^ Vt, samples information on a broader energy range 
around the Eermi energy, as pointed out recently—. 

To end, we have to admit that we are unable to unravel the 
origin of this mysterious ghost current that seems to be present 
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even in a standard vacuum nanogap. However speculative this 
idea might seem — a fact of which we are fully aware—, we 
believe that the experimental support presented here in favor 
of a presently unknown (“ghost”) transport channel can repre¬ 
sent a sufficient justification for considering the scenario pro¬ 
posed here a working hypothesis, which deserves further in¬ 
vestigations. Nanoelectronics could be the main benehciary. 
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Transition Voltage Simulation 




Fig. 3 (a) Transition voltage Vt and (b) Ohmic conductance G as a 

function of the nanogap size d for STM-tip of tungsten 

(<I>n' = 4.55 eV) and substrates of Pt at positive biases 

(<!>/), = 5.65 eV) in the presence of a ghost channel of dimensionless 

resistivity p = 10^. v-values are given in the legend. Image charge 

effects are included using the exact classical interactio n 

which is cutoff close to electrodes using the procedure described 

elsewhereii. 



Fig. 4 By subtracting larger and larger ghost currents -values 
at V = 1V specified in the legend) from the measured current, the 
measured (^-dependence of V) (black squares, same as the black 
circles in Fig. 4b of ref.—) switches from decreasing to increasing 
with 1 /d. The latter dependence is similar to that depicted in Fig.|^ 
expected without ghost contribution. (Courtesy of the authors of ref. 
0, who allowed to use here their experimental data.) 


Transition Voltage in OPI n Series 



Inverse molecular length 1/d (nm *) 


Fig. 5 The transition voltage of the OPI series exhibits a 
dependence on the molecular size (n)compatible to that found in 
Fig- S'- The points represent experimental results (cf. Table SI of 
ref. Ig). The error bars correspond to relative standard deviations in 
Vt of ~ 20%, as typically found in experiment s (Size ii = 5 has 
been eliminated since it seems to correspond to a special crossover 
situation—.) 
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